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Abstract 

We investigate the energy loss of quark and gluon jets in quark-gluon plasma produced in the 
central Au+Au collisions at RHIC energy. We use the physical characteristics of initial and mixed 



3- 



We also take into account the possibility 



phases found in the effective quasiparticle model 
of hot glue production at the first stage. We calculate the suppression of 7r° spectrum, which is 
caused by the energy loss of gluon and quark jets. We compare this suppression with data reported 
by PHENIX. We find that the suppression is described correctly by the quasiparticle model with 
decrease of thermal gluon mass and effective coupling in the region of phase transition (at T —* T c 
from above). Thus it is possible to investigate the structure of phase transition with the help of 
hard processes. We also show, that the energy loss at SPS energy is too small to be observable. 
This is caused, in fact, by too low initial temperature of the plasma phase at this energy. 
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I. INTRODUCTION 



Energy loss of high energy quark and gluon jets in relativistic A+A collisions leads to 
jet quenching and thus probes the quark-gluon plasma Recently the considerable sup- 
pression of hadron spectra at moderately high p± < 6 GeV was found in the central Au 
+ Au collisions, especially for it production [2j. This suppression can be explained by the 
energy loss of the moderate p± < 10 GeV gluons jets. This is one of the major new results 
at RHIC energy, which has not been observed before at lower SPS energy. The energy loss 
was investigated in the work [3] in various orders in opacity L/X g (where L is nuclear radius 
and X g is the gluons mean free path). It was shown that these series is strongly dominated 
by the first term. It was also shown in Ref. that finite kinematic boundaries decrease the 
energy loss at intermediate jet energies as compared to the asymptotic limit. Recently, in 
paper p], the physical characteristics of initial plasma phase and also of mixed phase at SPS 
and RHIC energies (i.e. the values T , r , r c ) were investigated on the basis of the quasipar- 
ticle model and isentropic expansion. In the present work we use these initial conditions for 
investigation of the energy loss of high energy gluon and quark jets at RHIC energy. In these 
studies we take into account also the possibility of two stages of equilibrium, i.e. production 
of the hot glue at first stage [6]. The calculation of energy loss is important for investigation 
of the phase transition structure of quark-gluon plasma into hadrons. In the Refs. [tIl! the 
phenomenological parametrization of coupling constant G(T) was used in accordance with 
new lattice data. When approaching to phase transition point from above, the decrease of 
thermodynamic values is caused, according to SU (3) gauge theory, by increase of the ther- 
mal gluon mass and also of the coupling constant G(T). In these conditions the applicability 
of perturbation theory is questionable. In phenomenological model of confinement jsj the 
decrease of thermal gluon mass m g (T) at T — > T c from above is assumed, which is connected 
with decrease of effective coupling strength G(T). It can be understood, as the more and 
more gluons become confined and form heavy glueballs with decrease of T, and the effective 
glueball exchange interaction between gluons are reduced. However, in this case the entropy 
density s(T) will exceed the lattice entropy because of domination of light masses near T c . 
This difference may be accounted for in quasiparticle model by modification of the number 
of effective degrees of freedom in thermodynamic functions: 

9 9 - C(T)g g . (1) 



The explicit value C(T) may be estimated as ratio of the lattice entropy and the quasiparticle 
entropy density s g (T) with decrease of the mass m g (T). At T 3> T c we have C(T) ~ 1, and 
near T ~ T c we have C(T) < 1. The value C(T) for gluons has the form P|: 

C(T,T C ) = C (1 + 6 C -^, (2) 

where Cq ~ 1.25, 5 C ~ 0.0026, (3 C ~ 0.31. The thermal mass m g (T) can be parameterized 
well by formula: 

N T 
rn g {T) = Q)^G T{l + 5-±Y, (3) 

where (3 ~ 0.1, 5 ~ 10~ 6 . It possible to vary some parameters, for example (3 ~ 0.05,5 ~ 



10 7 are reasonab 
from lattice data 



. The form (J3J)is analogous to behaviour of Debye mass ttld, extracted 
. The value Gq can be determined by asymptotic value of the thermal 



mass, chosen from coincidence of lattice and perturbative masses atT~3T c tals|. We had 



m 

found for pure gluons: G ~ 1.9 |5|. 

The relations ©, © with above-mentioned parameters give good description of SU(3) 
lattice data for entropy s, energy density e, and pressure p for T close to T c . It is possible to 
extend the effective quasiparticle model to a system with dynamic quarks, with analogous 
effective coupling G(T), and the function C(T,T C ) with some variation of parameters 0. 
One can show that the same constant Gn ~ 1.9 describes the lattice data for 2 and 2+1 
flavors (for example, the lattice pressure 0,0- 

It is interesting to note, that the massive constituent quarks (m q and m s ) appears also due 
to the decrease of number of degrees of freedom in the presence of octet of pseudogoldstone 
states. This is a consequence of conservation of the entropy and of the number of net 
nucleons. It can also be demonstrated that hadrons (in the hadronic part of the mixed 
phase) appears with the same effective number of degrees of freedom [5j. 

It can be shown, that such picture of phase transition (with decrease of effective coupling 
strength G(T) near T ~ T c ) gives reasonable quantitative description of the jet quenching 
and of suppression of the hadronic spectra in the central Au +Au collisions at RHIC energy. 
However, the use of the "running" coupling a s (T) in the perturbative decomposition of the 
thermodynamic values in plasma leads to too large energy loss of jets, which disagrees with 
experimental data for suppression of the hadrons with large p± at RHIC energy. 

In the Ref. P| we do not find a noticeable difference for meson and baryon spectra in 
the ordinary perturbative theory with the running coupling a s (T) and also in model with 



phenomenlogical parametrization of coupling G(T) (3, |8( in comparison with the effective 
quasiparticle model. However, the ordinary perturbative model disagrees with 577(3) lattice 
data in the region of phase transition. Thus the spectra of particles apparently weakly 
depend on the character of phase transition. Therefore the investigation of jet quenching 
represent significant interest. 

In Sec. Ullwe calculate the energy loss of the high energy gluon and quark jets in quark- 
gluon plasma at RHIC energy. In these calculations we use the initial conditions in plasma 
at SPS and RHIC energies, which were found in the quasiparticle model Q]. We also take 
into account here the production of hot glue at the first stage. We show that the energy loss 
at SPS energy is too small. We show also, that in the perturbative theory the energy loss 
of gluon jets is too large. 

In Sec. IIHI we calculate the suppression of tt° at moderate p± ( 3 < p± < 6) GeV in the 
central Au+Au collisions. We take into account the jet quenching and the parton shadowing 
factor in nucleus. We show that suppression of the 7r° spectra at PJ° > 3 GeV/c does not 
contradict the experimental data at RHIC energy. 



II. ENERGY LOSS OF HIGH ENERGY JETS IN QUARK-GLUON PLASMA 

In the effective quasiparticle model 5] we have investigated the initial condition and the 
evolution of the plasma stage, where there is equilibrium for both quarks and gluons. We 
have found for RHIC case the values T ~ 216.3 - 219.6 MeV, r ~ 2.22 - 2.18 fm and 
r c ~ 6.34 — 6.52 fm. We have also found corresponding values for SPS energy T ~ 175 
MeV, tq — 3.28 fm, r c ~ 4.1 fm, so, here we have the short plasma stage. The production 



re relatively large gg cross section 
In the lowest order, the matrix 



of a more hot glue plasma at the first stage is caused by t 
in comparison with the qg and qq cross section of Ref. [fj 
elements M 2 in formula: 

da -no 2 , , ^ 2 fA\ 

n * 

at large angles are related thus as [10]: 

K^JK^JK^ = 30.4/5.4/2.2 (5) 

i.e. gg scattering is most important. The small angle scattering leads to divergent cross- 
sections, which however are finite in QGP due to the finite "Debye mass" t m i n = g 2 T 2 . The 



cross section at large angles is defined by formula: 

30.4vra 2 95.46a? _ a 2 



large-angle ^ « ^ s ^ O (R\ 

The cross section at small angles has the form: 

small-angle ^ TO s _ ^f^£_ /y\ 

99 ~ 2 47rT 2 a s ~ 8 T 2 ' 1 J 

The effective scattering rate l/r 9 is determined by a sum of large and small angle values: 

I =ns (2.65^ + ^). (8) 

Here the value a s is determined by the effective coupling constant a s — f-, where: 

G(T,T c ) = G (l + 8-^y (9) 

and G ~ 1.9 (see Eq. ©) 

The net gluon density in effective quasiparticle model is: 

oo 

16T 3 f x 2 C(T,T c ) 

o e v 

where 

„n 2 



m 2 1 



67,(1 + 5-^ 



To estimate the initial temperature of the hot glue we use Bjorken model [ll[. The value 
dN/dy (number of charged and neutral particles per unit central rapidity) for Au + Au 
collisions can be described by formula 

(fL^ 8in ^ ;i - (i2) 

This value agrees with experimental value 1374 at \fs = 130 GeV. 

In this paper we make the following assumptions: at the first stage we have the hot glue 
in equilibrium and nonequilibrium quarks. The total entropy of gluons and such quarks (at 
temperature T g ) is less than the entropy of gluons and all quarks at more a later stage (at 
above mentioned temperature To). With expansion of the system the temperature decreases 
and the entropy of quarks and gluons decreases (if to assume that all the quarks are nonequi- 
librium ones). However, in reality, some part of quarks achieves equilibrium and entropy 



increases. With subsequent decrease of T the number of equilibrium quarks increases, and 
near To we have the equilibrium of both gluons and quarks. After that the cooling becomes 
isentropic. 

One can estimate the fraction of gluons and sea u, d, s quarks (and antiquarks) from 

the structure function 12], if to imagine the system moving backward in time from many 

secondaries into hot partonic plasma. This approach has to normalize parton multiplicity to 

the total entropy, later observed as multiplicity of secondaries. We assume its conservation 

at later stage (for T < T ) and its additional production at previous stage Q|. For average 

number of partons (for example, gluons in nucleon) we use formula [l^j ]: 

i 

N g = J dxg(x,Q), (13) 

m char /\p\ 

where g(x, Q) is structure function, m c har ~ rn p is characteristic hadrons mass, and \p\ is 

nucleon momentum. For RHIC energies we have m c h ar /\p\ = 2m c h ar / 'y/s rs 0.012. 

Similarly, the number of the sea quarks is 

i 

m char /\p\ 

This gives estimate N sea /N g « 0.44. This correspond to known fact, that gluons dominate 
in structure function at small x: the number of gluons is more than twice larger than the 
number of sea quarks and antiquarks. The characteristic average value x m \ n ~ 0.012 is 
close to value x m i n ~ 0.01 in the works of other authors (for example E.V.Shuryak). These 
meanings are reasonable for estimation of the hot glue temperature T g , since partons with 
x < x m i n carry only a small part of the total momentum. For estimation of the number of 
gluons in central region of rapidity we use the relation (fT2*|) : 



dN g dN sea dN va i 

— - H 1 

dy dy dy 

dN 9 , nAA dN 9 , dN ™l 



+ 0.44— + — ^ ~ 1374 , (15) 
dy dy dy 

where N va i = 3(N — N) ~ 49 (where N — N ~ 16.3 [5] is the number of net nucleons for 
RHIC in central region of rapidity). 



Hence, the estimate for gluons is: 



dN . 

" ] - 920. (16) 



dy 
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Taking into account Eq. (JHJ), we have the relation for determination of the hot glue initial 
temperature T g : 

dN q . x2 T 2 



n„r„ 



* 9 ' 9 dyitR\ u ~ 12 - 5 ^ m -) " 2.65a*(T) + 1.125a s (T) ' 

Here the physical values are in the units of m w = 139 MeV. The function a s (T) is determined 
by Eq. ©. We find from Eq. JTZJ): 

T g ~ 345MeV (18) 

From Eq. (JTUJ) the initial gluon density is: 

n g (T g ) ~ 19.42 m£. (19) 

Thus the time required to achieve the equilibrium for gluons is 

Tq ~ ~ — ~ 0.91 fm/c . (20) 

n g (T g ) m n 

The same time for quarks and gluons found in Ref. \^ is To — 1.54/m,,- ~ 2.18 fm/c. This 
agrees with estimate in Ref. p. 

In the region T < T < T g the value r g we estimate also by Eq. (JHJ). Such estimates are 
reasonable, as in the effective quasiparticle model the coupling constant a s (T) decreases with 
decrease of temperature, and Eqs. (JBJ)-(JHJ) in the lowest order are a good approximation. In 
fact, in the effective quasiparticle isentropic model we have the initial conditions in plasma: 
T ~ 219.6 MeV and r ~ 1.54/m^ [5]. However, from Eq. (JHJ) at T = 219.6 MeV and 
n g = 3.7ml (for the case of equilibrium quarks 0], see Table |TJ) we have n g r g = 5.67m 2 , 
which gives the same value r g = 1.53/m^. Thus Eq. (JHJ) provides matching of values r for 
the two mentioned regions. 

The values r for different values of T are shown in Table |U In this table we also give 
the estimates of equilibrium light quark density in the region T < T < T g , where the total 
equilibrium of gluons and quarks is not achieved (we assume that in the initial state of the 
hot glue at T g = 345 MeV all the quarks are nonequilibrium ones) 

These estimates can be done in the following way: 

1. We calculate the number of gluons N g = n g r g iiR\ u in the region of the hot glue 
T < T < T g . The number of the nonequilibrium quarks N q+S we estimate using Eq. (JTSj) 
and subtracting the value N g . 



2. For these quarks we calculate the nonequilibrium entropy S q + S in the Boltzmann ap- 
proach (the expected difference for calculations with fermions is about ~ 5%). We calculate 
also the entropy S g for equilibrium gluons. We assume, that the total entropy S g + S q+S 
increase smoothly (almost linearly like N q+S (T)) to the equilibrium entropy So(T ) ~ 6300 
with decrease of T — > T [5]. Using this assumption, we estimate the equilibrium addition 
from quarks to the total entropy density , and equilibrium addition n q (T) to quark den- 
sity in the region T < T < T g . The equilibrium densities n g {T) and n q (T) in Table |U (at 
T c < T < 7J) we find, using the formulas of effective quasiparticle model fl. 

The dominating first order radiation intensity distribution for expanding plasma is given 



in Ref . 



■Q: 



oo \k\i_ 

dl 9C R E 



J dz p(z) J <i 2 k , a & 



dx ir 2 

zo \k\ 

Qmax 



, 2 a s kq 

' d ,,2^12 1,2,, ~ X 





X 



1 — COS 



[q 2 + /i 2 (^)] 2 k 2 (k-q) 2 
(k-q) 2 (z-z ) 



(21) 



2x{l-x)E 

where E is jet energy, and Cr is color factor of jet (Cr = N c for gluons). Here we take into 
account substitution x — > x(l — x) for x — > 1 fl^ |. It is assumed that quark-gluon plasma can 
be modeled by well separated color-screened Yukawa potentials. The upper kinematic bound 



of medium induced momentum transfer is |g| m ax ~ a/3 / u(t)£'. The transverse momentum 
is connected with gluon emission distribution from parton jet in the absence of a medium. 
The kinematic bound on the transverse momentum is k 2 ^ = mm[AE 2 x 2 1 AE 2 x(l — x)} and 
kmin = A i2 ( r ) f° r gluons with light cone momentum fraction x (where /i 2 (r) = 47ra s T 2 (r)). 

The value z = t is limited by thickness L of target, but as a matter of fact, by end 
of plasma phase and by beginning of mixed phase r c , i.e. by duration of plasma phase 
t < r < t c . The value r c does not exceed the radius Ra u (or it is close to it). The value 
p(r) is the gluon or quark density at time r along the jet path (i.e. p = n g or n q ). The 



S 



energy loss AE is defined by integration of dl/dx in Eq. (|21|) over x. The value 

1 |k|max 

h(E,r(T)) = J dx J d 2 ka s (r) 

o |k| min 

| Q| max 

d 2 qa 2 (r) kq 
[q 2 + /i 2 (r)] 2 k 2 (k-q) 2 X 

(k - q) 2 r 

cos 



X 



(22) 



2x(l-x)E_ 

in the finite kinematic bounds is calculated by Monte Carlo method for values of t(T) in 
the region r g < r < r c (for RHIC energy) and for several values of E. In the calculations of 
complete energy loss AE the scattering of gluons (and quarks) by gluon and quark potentials 
in plasma medium should be taken into account for every value of E. For example, the 
fraction of gluon scattering on gluon potential in medium is defined by formula: 

""^ - ' ( 23 ) 



n g a gg + n q a gq 1 + ^ 

where a gg is the cross section for gluon scattered on gluon potential, and a gq is corre- 
spondingly cross section for gluon scattering on quark potential. Also the fraction of gluon 
scattering on quark potential is (1 + )~ 1 . Taking into account the decrease of the color 
factor we have u qg /a gg = 4/9 0. 

The complete energy loss AE of gluon jet is determined by formula: 

Tc 

where the value I q (E,t) corresponds to Ii(E,t) in variable \k\/fi and \q\/ fJ>. Similarly, the 
energy loss of quark jet is determined by formula: 



9n a 

T 9 



Tc 



l iThi E,r)^\ 
9ttV J uv ' ; i + ^a 



drIo(E,r)-^-, (25) 



's 



where Cf is the color factor of quark: Cp — ^ N . 

However it should be noted, that for the well-separated color screened Yukawa potentials 
the condition A = 3> 1/// must be implemented (here A is the mean free path of partons, 
and n is the color screening parton mass). This condition is well realized for parton scattering 
on gluon potential, where \j? ~ 47ra s T 2 , and we have A Si9 /i s ^> 1. But this condition is not 
realized for partons scattering on the quark Yukawa potential, where /i 2 = 4ira s T 2 /6, and 
we have here A ffi9 // ? < 1. Therefore the estimation of partons energy loss AE on gluon 
Yukawa potential is the most real in this model. These energy loss AE of gluon jet we 
give in Table |J for example at E = 10 GeV. But calculating also the energy loss on quark 
potential by formulas (j24l - !25jl . we find the decrease of AE no more than 10%. 

The complete energy loss AE is calculated by numeral integration over r. We have 
AE g ~ 3.15 GeV at E — 10 GeV. The complete energy loss of gluon jets AE g (E) is shown 
in Fig. ^opposite E (Fig. El shows the relative energy loss AE(E)/E). 
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Figure 1: The energy loss of gluon jet in quark- gluon plasma at RHIC energy. 
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Figure 2: The relative energy loss of gluon jet in quark-gluonf plasma. 

Let us consider now the energy loss at SPS energy. For the value dN/dy per unit central 

10 



Table I: The values of various physical quantities for expanding plasma at RHIC energy and energy 
loss of gluon jet AE. 



T (MeV) 


r (m 

V 7T / 


a s 


lo 

U 


n„ (rnt) 

y v 7r / 


n„ (m't) 


AE (GeV) 


345 


0.645 


0.267 


0.092 


19.42 





2.367 


325 


0.70 


0.265 


0.088 


15.95 


4.41 


2.097 


300 


0.79 


0.262 


0.085 


12.22 


4.71 


1.855 


275 


0.908 


0.26 


0.083 


9.09 


4.45 


1.585 


250 


1.075 


0.257 


0.081 


6.5 


4.21 


1.38 


225 


1.34 


0.252 


0.079 


4.4 


3.50 


1.13 


219.6 


1.54 


0.247 


0.078 


3.7 


6.1 


1.02 


210 


1.74 


0.241 


0.077 


3.12 


5.134 


0.964 


200 


2.04 


0.235 


0.075 


2.56 


4.195 


0.84 


190 


2.44 


0.229 


0.073 


2.045 


3.32 


0.778 


185 


2.72 


0.222 


0.067 


1.795 


2.9 


0.663 


180 


3.1 


0.215 


0.063 


1.545 


2.47 


0.589 


175 


3.7 


0.201 


0.055 


1.286 


2.03 


0.48 


172 


4.26 


0.184 


0.043 


1.116 


1.73 


0.37 


170 


4.6 


0.057 0.00067 


1.23 


1.61 


0.0213 



rapidity (charged + neutral) we have [5J]: 



- 803. (26) 

d V J y =o 

The number of partons from structure function is also obtained by Eqs. (JH1) , (HH) , but now 
x m i n ~ 0.09 — 0.1. The fraction of sea quarks is N sea /N g ~ 0.58. Taking into account the 
number of valent quarks in the central region of rapidity ~ 57 x 3 we have now: 



dN g 



~ 400. (27) 



dy 

We use the relation analogous to (fTTj) for derivation of the hot glue initial temperature T g : 

5.1. (28) 



1 fdN g 



TiR 2 Pb V dy 

It can be shown, however, that there is no solution of equation (fTTj) for T g in that case. 
Apparently, at such low energy the hot glue phase is practically indistinguishable from 
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Table II: The same as Table I but for SPS energy. 



T (MeV) r (m" 1 ) a s I n g {ml) n q (ml) AE (GeV) 



175 


2.4 


0.201 


0.061 


1.28 


2.2 


0.212 


174 


2.52 


0.197 


0.058 


1.224 


2,1 


0.204 


173 


2.65 


0.192 


0.055 


1.168 


1.99 


0.188 


172 


2.8 


0.184 


0.049 


1.11 


1.88 


0.168 


171 


2.9 


0.172 


0.04 


1.056 


1.76 


0.141 


170 


2.98 


0.057 0.00075 


1.23 


1.74 


0.0091 



phase of total equilibrium (or close to it). In Table HD the energy loss AE(t) is shown at SPS 
energy (E = 4 GeV) for the case of total equilibrium. At E = 4 GeV the total energy loss is 
AE ~ 106 MeV. Such small energy loss is caused by a short lifetime of the plasma phase at 
low SPS energy. It should be noted, that the energy losses, which are found here, correspond 
to the model of phase transition with decrease of effective coupling strength G(T) and mass 
m a (T) for T — > T c from above. 

n 

The similar problem was considered in the Ref. [5j using the ordinary perturbation theory 
for powers of running coupling a s in quark-gluon plasma (up to order 0(a s )). In this 
approach the coupling a s increases with T —>■ T c above. Although this model disagrees with 
SU(3) lattice data in the region of phase transition (i.e. close to T c ), we do not find here a 
noticeable difference for spectra of particles in comparison with the effective quasiparticle 
model (i.e. the spectra weakly depend on structure of phase transition). 

However, let us consider the energy loss at RHIC in this perturbative model. We use the 
Eq. (JHJ) with running coupling a s (with A ~ 180 MeV) and values n g without quarks and 
with quarks in the perturbative decomposition. We find the initial temperature of the hot 
glue to be T g ~ 400 MeV. The initial temperature of the total equilibrium is To ~ 219 MeV, 
t ~ 2.18 fm and r c ~ 3.73 fm 5]. The total energy loss can be found from corresponding 
Table, which analogous to Table |lj For example, we find AE ~ 13.5 GeV for E — 16 
GeV, AE ~ 8.5 GeV for E = 10 GeV, and AE ~ 3.5 GeV for E=4 GeV. These energy 
losses as too high and disagree with the data for suppression of hadrons with large p± in the 
central nuclear collisions at RHIC energy. This disagreement is connected with injustice of 

n 

perturbation theory in region of phase transition. In Ref. [5j| the quasiparticle model with 
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phenomenological parametrization of coupling constant G(T) was also considered. In this 
model G(T) increase for T — > T c from above. This model agrees with the new lattice data 
and provide also good description of the baryon and meson spectra. However, it can be 
shown that in this model we have a paradox result for energy loss of gluon jet in plasma: 
the energy loss exceeds the energy of jet itself (even without taking into account the effect 
of the hot glue). This contradiction is caused by too large value of coupling G(T), especially 
close to the phase transition point. 



III. SUPPRESSION OF PIONS WITH LARGE p ± IN CENTRAL AU+AU COL- 
LISIONS 

The jet quenching reduced the jet energy before fragmentation, where the jet transverse 
momentum is shifted by energy loss on to value AE(E) 15]. To account for this effect we 
should replace the vacuum fragmentation function by effective one z*/z c D h / c (z*, Q 2 ), where 



1 E 

The invariant cross section of hadron production in central A+A collisions is given by: 

b max 

r AA 



drr AA I" 
E h^jjh= J d 2 bd 2 rt A {r)t A {\h-v\) 



o 

J dx a dxbd 2 k± :a d 2 k± t i J x 

abed 

x 9A{k±, a , Q 2 , r)g A {k±,b, Q 2 , |b - r|) x 
x fa/A(x a ,Q 2 ,r)f b/A (x b ,Q 2 ,\b - r|) x 
x daz* D h/c (z* c ,Q 2 ) 
di z c irz c 

Here t A (r) is the nuclear thickness function, kx,a an d are the initial transverse mo- 
menta of partons, f a / A and f b /A are parton structure functions. It is usually assumed that 
distribution g A (k±) has a Gaussian form. 

It should be noted, that intrinsic k± and the transverse momentum broadening (Cronin 
effect) are important for final hadron spectra for SPS energies. However, with the increase 
of energy the spectra become flatter and small initial k± correspond to small variation of 
spectra. At RHIC energy one can neglect the effects of initial k±_ with a good precision 
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16| . In this work we do not take into account the intrinsic transverse momentum at RHIC 
energy. 

For parton distribution of nucleon in nucleus we take into account the parton shadowing 
factor S a ^/A(x, r) for spherical nucleus and S a ^/A{x) for flat disk, for which we take the 
parametrization used in HIJING model [l2|. In this model the factor S a /^(x,r) is splitted 



into two parts: S(x,r) = Sq(x) — ctA(r)Si(x), where ^(r) = 0.1(A 1 ^ 3 — 1)4/3a/1 — r 2 /R 2 A . 
For flat disk we have ^(r) — > a a = 0.1 (A 1 / 3 — 1). 

The upper limit for the impact parameter is 6 max ~ 0.632 Ra u for 10% central Au+Au 
collisions. Neglecting initial k± in Eq. ()30j) it is convenient to introduce the new variable 
x a = X\ = x' n L £/2z, x b = x 2 = x^£/2z(£— 1), where x\ = 2p K L / y /s, z = x\jx^ x± = 2E/y/s, 
i.e. z = p±/E. After integration over r and b in Eq. ()30)) we obtain cross section for hadron 
production in gluon jet in the central Au+Au collision: 



drr AA QK f „ 1 

Eh ^i = My J dz^(Q(z))z^D h/g (z:,Q 2 ) {1 _^ {z) - 



2z 

x ± 



x J dtf™(Z)<!>(x 1 ,x 2 )x 1 G g (x 1 )x 2 G g (x 2 ), (31) 

2z 

where Q = p±/2z, Q = p±/2z*, and G g (xi), G g (x 2 ) are the gluon structure functions. Dh/ g 
is fragmentation function 18| and K is K-factor. Function <&(xi,x 2 ) describes the contribu- 
tion of shadowing. For spherical nuclei (for 10% central collisions) we obtain: $>(xi,x 2 ) = 
226.58 3 (xi)S Q (x 2 ) - 117.27 S^xJSoip) - 97.69 S (x 1 )S 1 {x 2 ) + 54.87 S^xJS^). Func- 
tions 5*o and Si can be found in Ref. [17]. The value f gg (£) is the fraction of gluon-gluon 
elementary cross-section in variable £: 

f (e) gce^i) (e-i) 2 i (e-i) 2 m 



We also have: 



f (c) ! + ^ 2 1 4 (e-i) 1 4(e - 1) , 33) 

= l^pry ( 34 ) 
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The value z mSLX can be found from Eq. (|2*5j) : z* = Px Afi = E V _\ E = 1 (taking into ac- 

\ E ' 

count that z = p\jE\ The value of E which corresponds to gluon jet with some p\ can 
be found from Fig. ^in large scale). For example, for values p\ = 3, 4, 6 GeV/c we can 
estimate z max = 0.61, 0.64, 0.67 (e.g. at p ± = 4 we have E - AE(E) ~ 6.25 - 2.25 = 4, and 
z max ~ 4/6.25 = 0.64). 



We obtain the approximate numeral formulas for ^§-{z) in the region z min < z < z mSLX 
for several values of p\ (which are (3 — 5)% precise). For example, for p\ = 4 GeV we have 
(here z min = 2Pjy/s ~ 0.0615): 

AE 

\z) ~ 0.46z + 0.061 + 0.95(z - 0.0615)(0.64 - z). (35) 



The relations of this kind we use in formula ()31|) for gluon jets. 

But it should be noted, that for the quark scattering on gluon and quark potentials 
the color factor decrease 16/81 = 0.198 times and the energy loss AE should decrease 
correspondingly. From relation E — O.lOSAi*^ = p\ we can find for p\ = 3, 4, 6 GeV the 
corresponding values z max ~ 0.91, 0.92, 0.924. Here we also find the approximate analytical 
formulas for ^jjr{z), which we use in Eq. ()31j) for quark jets. For example at p\ = 4 GeV we 
have: 

-=-(z) ^ O.O682 + 0.014 + 0.11(2 - 0.0615)(0.92 - z) (36) 
E 

In this work we calculate the effective suppression factor for the 7r° spectra at RHIC 
energy, or the ratio: 

dN AA 



*-<pj = • < 37) 

^in dyd 2 p ± 

frmax / ^max 

where T = J T(b)d 2 b / J d 2 b. This is ratio between the spectrum in central AA and 
o / o 

pp collisions, which is normalized to the effective total number of binary NN collision in 
central AA collisions. If we do not take into account the nuclear effects (jet quenching and 
shadowing), this ratio should be unity at large p±. 

In Eq. (13 lj) for invariant cross section we take into account the following parton collisions: 



1. The gluon-gluon collision, the gluon scattering by gluon potential and gluon fragmen- 
tation into 7T°. This process is described by Eq. (|3ip with effective fragmentation 
function D^o/g^z*, Q 2 ) and with function ^r(z) of the kind ()35|) for each value of 

Pi- 
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2. The gluon-quark collision (with u and d quarks), the gluon scattering by gluon and 
quark potentials and fragmentation of gluon. In Eq. (}3*T]) we use f qg (£) from (J3l?)) with 
corresponding factor, the product G g G q from and AE/E of the kind 

3. The qg collision, the quark scattering by gluon and quark potentials and fragmentation 
of quark into it . In Eq. (|3ip we have now D^a/q jl8| with ^r(z) of the kind (}3l)|) for 
each value of p\ due to the decrease of the color factor. 

4. The quark-quark collision (jHI|), the quark scattering by gluon and quark poten- 
tials and fragmentation of quark. We have here AE/E also of type (pjj) and product 

G q G q . 

5. The collision of gluon with sea quarks, the gluon scattering by gluon and quark poten- 
tials and gluon fragmentation. Here we also use the effective fragmentation function 
DjM/g with AE/E of the kind (|35j) and corresponding structure function from Ref. 

Q 

In order to find the effective suppression factor Raa(p±) of the 7r° spectra, the cross 
section (J31)) should be divided by the same cross section for the same value of p'f , but 
without nuclear effects (AE = 0,z max = 1), and without shadowing (Sq(xi) = S (x 2 ) = 1, 
Si(xx) = Si (xz) = for spherical nuclei). 

We have calculated such ratio for hard collisions (3 GeV/c < p'f < 6 GeV/c). For 
example, this ratio for p 7 ^ = 4 GeV/c (the contributions of collisions listed above in items 
1-5 are given in numerator and denominator) is: 

Raa(pH = 4GeV/c) = 

0.37 + 0.08 + 1.09 + 0.27 + 0.015 n ^ /oc , 

= ~ 0.27. 38) 

4.0 + 0.53 + 1.8 + 0.40 + 0.161 v ; 

The values of Raa corresponding to p^f = 3, 4, 5, 6 GeV are Raa = 0.23, 0.27, 0.32, 0.38. 
We show in Fig. El the ratio Rau,Au(p±) and the experimental data at RHIC energy. The 
theoretical values of Rau,Au are somewhat smaller, than the experimental values, though 
they apparently are within the limits of experimental uncertainly. 

Theoretical values of Raa in Fig. 01 correspond to sufficiently large values of p±. In the 
intermediate region 1 < p± < 3 GeV there are uncertainties associated with the interplay of 
contributions from hard and soft processes. In the soft region p± < 1 GeV/c it is possible 
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Figure 3: The suppression factor Raa(p±) for charged hadrons and neutral pions. The gray band 
for 7T° is defined by the sums of squares of the systematic errors of the measurement and uncertainty 
in the N + N reference and in {Nbi nary ). The errors for (h + + h~)/2 are not shown. The line 1 is 
the calculated factor Raa{p^) for hard processes, the line 2 is factor for soft processes calculated 
in assumption that the low p± spectra scale like A . 

make the estimate of R A a by Eq. (JHTJ) if the scaling A 1 - 1 is used: ~ ^l.i^g^ This 

gives the value Raa — 0.25 — 0.29, which is close to experimental data available for charged 
particles. 



IV. CONCLUSION 



In this paper we investigate the energy loss of quark and gluon jets in the quark-gluon 
plasma and suppression of the ir° spectra at RHIC energy. We use the initial conditions 
in plasma, which were found earlier in the effective quasiparticle model for SPS and RHIC 
energies p. We take into account also the possibility of production of the hot glue at the first 
stage with more high temperature and density. It can be shown, that plasma is sufficiently 

Tc 

thin: n = J dr n g (T)a gg (T) ~ 1.28, where n is an average number of jet scatterings. The 

energy loss in expanding plasma is calculated in the dominant first order P, |J] taking into 
account the finite kinematic limits. The energy losses AE of gluon and quark jets are 
calculated in a wide range of parton energies E. These energy losses are used for calculation 
of suppression of 7r° spectra with moderately high p± (3 < p± < 6 GeV/c), which is done 
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by modification of gluon and quark fragmentation functions. We also take into account 
the parton shadowing factors in nucleus 17^ . We obtain a considerable suppression of tt° 
spectra, which is caused by the effects mentioned above. This suppression agrees with the 
data reported by PHENIX in the region = 3 — 4 GeV/c. The estimates of Raa (I37|) 
in the soft region (p± < 1 GeV/c) also agrees with experimental data if the scaling A 11 is 
taken into account. 

The intrinsic transverse momentum k± at RHIC energy is neglected in this work. It is 
important at SPS energy (Cronin effect), but at RHIC energy one can neglect the effects 
of initial k± with a good precision Q|. In principle a small increase of Raa is possible at 
moderate p± because of transverse momentum broadening. We are planning to take into 
account the intrinsic k± in the next work. 

The most important conclusion of this work is the possibility to investigate the structure 
of the plasma phase transition into hadrons with the help of hard processes. The correct 
quantitative description of 7r° suppression is probably possible only in a model of phase 
transition which includes the decrease of thermal gluon mass m g (T) and effective coupling 
G(T) near the phase transition (T ~ T c ) point. However, the quasiparticle model with 
increase of these values at T —>■ T c leads to a too high energy losses (which can even exceed 
the energy of the jet itself). The main reason of this problem is too large value of the 
coupling G(T) near T c . We show also that energy losses at SPS energy are very small as a 
consequence of a too low initial temperature of the plasma phase at this energy. 

It is interesting to investigate also the suppression of charged meson and baryon spectra, 
which we are planning to do in the forthcoming paper. 
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